Spin chirality on a two-dimensional frustrated lattice 
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The collective behavior of interacting magnetic moments can be strongly influenced by the topol- 
ogy of the underlying lattice. In geometrically frustrated spin systems, interesting chiral correlations 
may develop that are related to the spin arrangement on triangular plaquettes. We report a study of 
the spin chirality on a two-dimensional geometrically frustrated lattice. Our new chemical synthesis 
methods allow us to produce large single crystal samples of KFe3(OH)s(S04)2, an ideal Kagome 
lattice antiferromagnet. Combined thermodynamic and neutron scattering measurements reveal 
that the phase transition to the ordered ground-state is unusual. At low temperatures, application 
of a magnetic field induces a transition between states with different non-trivial spin-textures. 



Geometrically frustrated magnets are unusual in that 
they may have disordered ground-states in which an 
enormous number of spin configurations share the same 
energy. 0,0 The Kagome lattice antiferromagnet, formed 
of corner sharing triangles, is one of the most highly 
frustrated two-dimensional (2D) lattices. Though long 
regarded as a prime model for studying spin frustra- 
tion, Q, lUElOJ |o] the Kagome lattice has escaped precise 
magnetic characterization because compounds compris- 
ing this lattice are difficult to make pure and in large sin- 
gle crystal form. We present a study of an ideal Kagome 
lattice compound, the iron jarosite KFe3(OH)g(S04)2. 
This material possesses robust chiral correlations related 
to the arrangement of spins around triangular plaque- 
ttes. [21 12( The presence of spin chirality in condensed 
matter systems may play a role in important phenomena 
ranging from high-temperature (high -TV ? ) su perconduc- 
tivity |9( to the anomalous Hall effect [lCllllj . Currently, 
there are relatively few experimental studies of spin chi- 
rality in frustrated magnets. Here, we study the non- 
trivial spin-textures related to chirality in a jarosite ma- 
terial in both the spin-ordered and spin-disordered states. 
Our magnetic susceptibility, specific heat, and magnetic 
neutron scattering measurements on pure single crystals 
show that the transition to the low temperature ordered 
phase is unusual. In addition, we find that the spin- 
texture of the ground state can be controlled by applying 
a magnetic field. 

The magnetic Fe 3+ ions of jarosite have spin-5/2 and 
lie at the corners of the triangles of the Kagome lat- 
tice, shown in Figure 1. To first approximation, the spin 
Hamiltonian is given by 
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FIG. 1: The Kagome lattice with spins arranged in two dif- 
ferent configurations, a, The "q=0" structure, which is the 
ground state configuration for iron jarosite. The spin arrange- 
ment has uniform, positive vector chirality, indicated by the 
+ within each triangular plaquette. b, An alternate spin 
arrangement with staggered vector chirality, known as the 
structure. 



where J is the nearest neighbor magnetic exchange in- 
teraction (J > 0). Because the interaction is antiferro- 
magnetic, the spin system is frustrated, and the corner- 
sharing arrangement leads to a higher degree of degen- 
eracy for the ground state than the triangular lattice. 
The iron jarosite is a particularly ideal Kagome lattice 
compound for the following reasons. First, it consists 
of single layers of undistorted Kagome planes, and these 
planes remain undistorted down to the base temperatures 
of our measurements (T~5K). Second, this jarosite can 
be synthesized with compositions that are stoichiometri- 
cally pure, as we describe below. This ensures that we 
are primarily studying the effects of geometric frustration 
rather than the effects of disorder. Third, large single 
crystals can be made, which allow investigation of the 
spin correlations of this Kagome compound that would 



2 



not be possible with powder samples alone. 

Until recently, jarosites have been prepared typically 
by precipitation under hydrothermal conditions: 

3Fe 3+ + 2K 2 S0 4 + 6H 2 0^ 

KFe 3 (OH) 6 (S0 4 )2 + 3 K+ + 6 H+. (1) 

Under these conditions, the monovalent K + cation is sus- 
ceptible to replacement by hydronium ions and the cover- 
age of the Fe 3+ lattice sites is incomplete. Also, only mi- 
crocrystalline materials are obtained owing to the height- 
ened acidity of the solution as well as the sjjeed and in- 
tractability of the precipitation reaction. E] The chal- 
lenges confronting the synthesis of pure jarosites have 
been overcome with the development of redox-based hy- 
drothermal methods. Control over the precipitation 
of the jarosite is achieved by inserting two oxidation- 
reduction steps before reaction (1) 

Fe + 2 H+ -> Fe 2+ + H 2 (2) 
2Fe 2+ + i 2 + 2H+^2Fe 3+ + H 2 0. (3) 

In this manner, the Fe 3+ is slowly generated through- 
out the course of the hydrothermal process and the pH 
is moderated because three overall equivalents of pro- 
tons are consumed in the production of an equivalent of 
jarosite. We have refined the synthesis process to opti- 
mize the size of the single crystals. Crystals as large as 
10 mm in length and 48 mg in mass have been grown, 
making possible the inelastic neutron scattering measure- 
ments that we discuss below. Details of the single crystal 
synthesis may be found in the Supplementary Informa- 
tion. 

Even though the Kagome lattice antiferromagnet 
should not order at any non-zero temperature, powder 
neutron diffraction measurements on KFe3(OH)6 (804)2 
indicate that the spins order in a coplanar 'q=0' arrange- 
ment below Tjy ^ 65 K. 0,0 In the q=0 structure, the 
spins on each triangle are oriented at 120° to each other, 
and the 2D magnetic unit cell is identical to the 2D struc- 
ture unit cell, as shown in Fig. la. The ordered spins can 
be decomposed into three sublattices, with the spins on 
each triangle labeled as S±, S2, and S3. The vector chi- 
rality for each triangle may be defined as: 

K v = ^=0& xS 2 + S 2 xS 3 + S 3 x Si). (4) 

For the coplanar arrangement, this vector is parallel to 
the c-axis with amplitude +1 or -1. The neutron pow- 
der results 01 indicate that each triangle has positive 
chirality (+1) in the ordered state, such that the spins 
point directly toward or away from the center of each 
triangle. 

For this compound, like as for real magnetic materials, 
the interaction Hamiltonian of the spins contains terms 
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FIG. 2: Magnetization and specific heat measurements of 
KFe3(OH)6(S04)2- a, M/H versus temperature with the ap- 
plied field in two orientations, H || c and H || ab measured on 
a single crystal sample using a SQUID magnetometer. The 
inset shows M/H for H \\ c for fields up to 14 Tesla mea- 
sured using an AC coil set. The upturn in the magnetiza- 
tion at low temperatures (below 20 K) can be described by 
Curie behavior in which the density of free spins corresponds 
to ~ 0.3% of the total number of spins, b, Specific heat of 
powder samples of KFe3(OH)6(S04)2 (closed symbols) and 
the non-magnetic isostructural compound KGa3(OH)6(S04)2 
(open symbols). The KGa3(OH)6(S04)2 data are scaled to 
match the Fe-jarosite data at high temperatures and are used 
to estimate the phonon contribution. 



beyond isotropic Heisenberg exchange. These additional 
terms cause the system to order at a non-zero Tjv and 
determine the ground-state spin arrangement. Currently, 
it is not clear whether the transition to long-range order 
is driven by weak intcrplanar coupling, spin-anisotropy, 
anisotropic exchange or some combination of these. An 
ordered phase on a Kagome lattice is characterized by 
two order parameters (the sublattice magnetization and 
the vector chirality) which have different symmetries. An 
intriguing possibility is that these symmetries are broken 
at different temperatures. This has been proposed as 
a result of numerical work on XY or planar triangular 
lattice systems. 0, However, this has not been con- 
clusively observed by experiment, (ltl Il7| For Kagome 
lattice systems, relatively little is known about the na- 
ture of the phase transitions that occur. 
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Magnetization and field-induced transition 

To address the above issues, we have first measured 
the magnetization of a single crystal sample (mass of 13.5 
mg) of KFe 3 (OH) 6 (S04)2 using a SQUID magnetometer. 
Measurements were taken with the applied field oriented 
along the c-axis (H || c) and within the ab plane (H \\ ab) 
as shown in Fig. 2a. At low fields (H < 5 T) with the 
field along c, a sharp peak appears near 65 K, indicative 
of the transition to the 3D magnetically ordered state. 
This is consistent with previous measurements on powder 
samples prepared under similar synthesis conditions. (Isf 
When the field is aligned along the afe-direction, the sharp 
peak is absent and is replaced by a broad cusp. 

In Figure 2b, we plot the specific heat C of powder 
samples of KFe3(OH)6(S04)2- A peak in the specific heat 
is found at the magnetic transition temperature of 65 K. 
The entropy associated with the 3D magnetic transition 
(integrating C/T over the temperature range from 2 K 
to 100 K) represents ~ 50% of the Rln6 (where R is the 
molar gas constant) total entropy expected for the spin 
5/2 system. This suggests that short-range correlations 
have already formed at much higher temperatures. 

At high temperatures (T > 150 K), the susceptibility is 
isotropic and follows a Curie- Weiss law x = C/(T—6aw) : 
consistent with previous results on powder samples. |l8| 
Fits to this law between 150 K and 550 K yield the val- 
ues C w = -800(30) K, and C = 5.6(2) cm 3 K/mol 
Fc. Because the data are taken for T < \9cw\, we ex- 
tract the effective moment /i e g and the nearest neighbor 
exchange coupling J using the high-temperature series 
analysis of Harris et.al. for the Kagome lattice. Q Our 
results indicate that J = 45(2) K = 3.9(2) meV and 
fi e g = 6.3(2) /ib (close to the spin-only value of 5.92 \xb 
for Fe 3+ ). We note that a small next-nearest neighbor 
interaction (J 2 > 0) would serve to reduce the calculated 
value of J by an amount of order 

The peak in M/H for H \\ c at T = 65 K indicates the 
presence of weak ferromagnetism along the c-direction. 
Recent theoretical work shows that antisymmetric ex- 
change, via the Dzyaloshinsky-Moriya (DM) interaction, 
may induce such a moment by canting the spins slightly 
out of the plane. This interaction is present if there is 
no inversion center between magnetic ions and adds the 
term 

E 3 ■ x Si) 

(id) 

to the spin Hamiltonian, where D is the DM vector. [2(| 
This interaction causes the spins on each triangle in the 
jarosite to form an 'umbrella' structure and gives each 
Kagome plane a net ferromagnetic (FM) moment. How- 
ever, the intcrlayer coupling in iron jarosite causes the 
FM moments to couple antiferromagnetically between 
layers in the absence of an applied field, as shown in 



the inset of Fig. 3b. Our measurements on single crys- 
tals allow us to explore this model of canted moments. 
The inset of Fig. 2a shows M/H as a function of temper- 
ature measured in high fields (H > 7 T). At these high 
fields, the peak broadens and the downturn in the mag- 
netization shifts to lower temperatures. Such behavior 
has been observed in the square-lattice antiferromagnet 
La2Cu04 |2i| for which the DM interaction does indeed 
result in weak ferromagnetism. [22^ 

To investigate this spin canting further, we made mag- 
netization measurements as a function of magnetic field 
along c direction as shown in Fig. 3a. The results show 
an abrupt change in the magnetization at a critical field, 
He, which we define as the field at which dM/dH is a 
maximum. We interpret this abrupt increase as a change 
from canted moments being oppositely directed between 
planes to canted moments aligned in the same direction. 
This is most likely caused by a 180° rotation of all spins 
on the layers that were previously oppositely canted, as 
shown in the inset of Fig. 3b. The critical field as a func- 
tion of temperature is also shown in Fig. 3b. For com- 
parison, we plot the integrated intensity of the (1,1,1) 
magnetic Bragg peak measured with neutron diffraction 
in Fig. 3c. We find that He scales quite closely with 
the staggered moment (which is proportional to the 
square-root of the Bragg intensity). The mean-field re- 
sult for La2Cu04 gives H n oc M< fx 1 where x^ is the 
2D staggered susceptibility. [2J| In our case, the staggered 
susceptibility for the Kagome lattice is expected a have 
weaker temperature-dependence than the square lattice 
due to the geometrical frustration. j2^| 

Below T — 49 K, the critical field becomes larger than 
14 T (the maximum field of our magnetometer). Hence, 
it is difficult to extract the zero-temperature values for 
He and the canted moment. However, an estimate for 
the canted moment can be made by considering the jump 
in the magnetization in the vicinity of He- For the data 
taken at T = 50 K, we obtain a canting angle for the 
ordered moment of 0.65(6)° with respect to the Kagome 
plane for H ~ Hq (this angle will increase upon cooling, 
as the order parameter has not yet reached its low tem- 
perature value) . The field-induced transition results from 
a competition between the interlayer coupling J± and the 
Zeeman energy; the magnitude of J± may be estimated 
from the following relation: fl" o (0)Mp(0) = 2S 2 \J±\, 
where Hc(0) and Mp(0) are the critical field and FM 
moment per Fe atom, respectively, at T = 0. Our re- 
sults indicate |Jj_| = 0.007(6) meV where the large error 
bar comes from the uncertainty in extrapolating He and 
Mf to T = 0. We find that the magnitude of J± is sev- 
eral hundred times smaller than the nearest-neighbor J, 
attesting to the two-dimensionality of the system. 

The field-induced spin-canting transition corresponds 
to a non-trivial change in the spin-texture of the jarosite 
sample. In particular, the transition yields a net, non- 
zero value for the scalar chirality, defined on each trian- 
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FIG. 3: Measurements of the field-induced transition to a 
state with non-zero scalar chirality . a, Magnetization versus 
applied field for several temperatures with H || c. The He 
values indicate the fields where dM/dH is maximum, b, He 
versus temperature. The inset depicts the change in the spin 
configuration below and above He- c, Integrated intensity 
of the (1,1,1.5) magnetic Bragg peak measured with neutron 
diffraction on a single crystal in zero field. This quantity is 
proportional to (M^) 2 , the square of the staggered moment. 
The lines in panels B and C correspond to power-law forms 
for H c and M 1 " with the exponent (3 = 0.25 and T N = 64.7. 



gular plaquette as 



K s = Si ■ (S 2 x S 3 ). 



(5) 



The presence of this type of chirality (in static or fluctu- 
ating forms) can have important consequences in strongly 
correlated electron systems, such as yielding an anoma- 
lous Hall effect in metallic materials. [l0j,|ll| For H < He 
the net scalar chirality for our jarosite sample is zero 
because the contributions from neighboring planes are 
equal and opposite. However, for H > He the spins flip 
(as depicted in Fig. 3b), and the net scalar chirality be- 
comes non-zero. There are few materials with non-zero 



scalar chirality in the ordered state, especially on a two- 
dimensional lattice. In iron jarosite, we have discovered 
a phase transition in which a net scalar chirality can be 
'switched on' by a magnetic field. 

Neutron scattering and spin fluctuations 

To probe the microscopic behavior of the magnetism 
in the temperature regime above Tjy, we made inelastic 
neutron scattering measurements of the spin fluctuations 
at T = 70 K within the L = plane of reciprocal space. 
For the temperature range T/v <T < 120 K, the suscep- 
tibilities for H || c and H \\ ab deviate from each other, 
indicative of the growing influence of spin-anisotropy. For 
the 3D ordered spin structure below Tjv, the stacking ar- 
rangement of the planes doubles the magnetic unit cell 
with respect to the structural unit cell such that the 
magnetic Bragg peaks occur at half-odd-integer values 
of L (L^0). However, for temperatures above Tjv, the 
correlations between layers are destroyed, and the 2D 
spin fluctuations yield "rods" of scattering in reciprocal 
space along the L direction. We have verified that all 
of the signal in Fig. 4a disappears upon cooling below 
Tjv, as expected for the transfer of the intensity of the 
2D critical scattering into 3D Bragg points lying out of 
the scattering plane. Therefore, the scans in Fig. 4a pass 
through these 2D rods of scattering and directly measure 
the dynamic structure factor of the spin correlations of 
the single Kagome planes in iron jarosite. The instanta- 
neous spin correlation length, measured separately in an 
energy-integrating configuration, is £ = 20(2) A at this 
temperature. A reciprocal space map of the intensity of 
the spin fluctuations is shown in Figure 4a. Outlined in 
red are the boundaries of the structural Brillioun zones 
(BZ). The strongest scattering occurs at the centers of 
certain BZ. The neutron scattering intensity can be cal- 
culated by 



I cx f(Qf 



QaQp) S al3 (Q,uj) (6) 



where a, [3 refer to x, y, z vector components, f(Q) is the 
magnetic form factor, and S(Q, ui) is the space- and time- 
Fourier transform of the spin-spin correlation function. 
The intensity variation reveals a wealth of information 
about the short-range ordered state in the critical regime 
just above T/y. 

First, the fact that the intensities are centered within 
the structural Brillioun zones indicates that the fluctu- 
ations have the q=0 arrangement. The absence of scat- 
tering at the (2,0,0) position is consistent with this spin 
arrangement. Second, the nearly equivalent intensities 
at the (1,0,0) and (1,1,0) positions reveal that the spin 
direction on each sublattice is not fixed. Rather, the 
two orthogonal spin directions within the Kagome plane 
are equally likely. Hence the spins fluctuations have XY 
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FIG. 4: Inelastic neutron scattering data for KFe3(OH)6(SC>4)2 measured above Tat, along with structure factor calculations. 

a, Intensity contour plot of data from inelastic neutron scattering measurements of a single crystal sample (mass=48 mg) . The 
bright regions are the loci of scattering intensity for the low energy [to — 1.5 meV) spin fluctuations at T — 70 K above Tn- 

b, Model calculations of the intensity as described in the text. The left plot corresponds to short-ranged q — correlations 
with uniform vector chirality, whereas the plot on the right depicts \/3 x y/3 correlations. Both calculations are for the case in 
which there is no preferred spin direction within the Kagome plane. 



symmetry at this temperature (the spin-rotational sym- 
metry within the Kagome plane is not broken). In con- 
trast, at low temperatures in the 3D ordered state, a 
preferred spin direction is chosen, and the intensities at 
(1,0,L) and (1,1,L) differ considerably. [HH 

The observation of short-range q=0 correlations in the 
2D fluctuations implies a particular arrangement of the 
vector chirality. In each region of correlated spins, the 
vector chirality must be uniform (all positive or all neg- 
ative for each plaquette). Using equation (6), we calcu- 
lated the neutron scattering intensity arising from q=0 
correlations between coplanar spins assuming only pos- 
itive chirality. Our model calculation are based on 7- 
unit cell clusters of spins, and we have averaged over all 
spin directions within the Kagome plane. The results are 
shown in Figure 4b, and the agreement with the data is 
excellent. Also shown in Fig. 4b is a calculation for the 
intensity if the chirality were staggered as in the v^x V3 
arrangement. If such correlated regions exist, the error 
bars on our data indicate that the fraction must be less 
than 5%. 

Our results shed light on several basic questions re- 
garding the magnetic phase transition in iron jarosite. 
First, we find that the instantaneous spin correlations 
above Tn are two-dimensional in nature, and the q—0 
arrangement is preferred. Therefore, the selection of g=0 
order (as opposed to \fi x y/3 order which is predicted 
to be the preferred ground state for the pure Heisenberg 



model |4i, 25] ) is caused by interactions within a single 
Kagome layer and is not controlled by the interplanar 
interaction |26| . Second, our neutron measurements re- 
veal critical spin fluctuations above Tjv that have XY 
symmetry; hence^the magnetic ordering is not driven by 
2D Ising physics [3 • Most interestingly, we find that the 
spin-rotational symmetry and the vector chiral symmetry 
are not broken simultaneously at T/y. 

The presence of vector chiral order above Tjv may be 
naturally explained in light of the DM interaction. One 
possibility is that the DM interaction is the dominant 
source of spin anisotropy in iron jarosite. In this case, 
the vector chiral order appears concomitantly with the 
growing spin correlations and does not represent a spon- 
taneously broken symmetry. Another possibility is that 
the XY anisotropy has an origin (such as symmetric ex- 
change anisotropy [2]} ) distinct from the DM interaction. 
Once the spin correlations become coplanar within the 
Kagome plane, the vector chirality of a triangular pla- 
quette becomes a discrete symmetry (the chirality vector 
is either up or down with respect to the c-axis). Hence, 
long-range chiral order is not precluded by the Mermin- 
Wagner theorem. [2^ However, chiral order is easily dis- 
rupted on a Kagome lattice by the proliferation of do- 
main walls (thermally induced defects) which can form 
with little cost in energy. 0, |2^, |2^| Then the presence 
of a small DM interaction (in particular, the non-zero 
out-of-plane component D z ) selects a particular chiral- 
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ity for each triangle, and thereby inhibits domain wall 
formation. It remains possible that the vector chirality 
represents a spontaneously broken symmetry; however, 
clarification of this point requires further neutron scat- 
tering measurements at higher temperatures. 

Our measurements of this ideal Kagome material re- 
veal new magnetic behavior related to two types of spin 
chirality: vector and scalar. The DM interaction is a 
significant perturbation to the Heisenberg Hamiltonian 
and strongly influences the low temperature physics. For 
T > Tjv, the vector spin chirality is ordered even in the 
absence of broken spin-rotational symmetry. In the or- 
dered state below T/v, we have discovered a field- induced 
transition to a state with non-zero scalar chirality. Thus, 
materials based on jarosites may be promising candidates 
for studies of the coupling between non-trivial spin tex- 
tures and the transport of electrons in frustrated systems. 
For example, carrier-doped compounds would probably 
show an anomalous Hall effect of topological origin, |HJ 
and this might have useful applications in spin-based 
electronics. 

Methods 

The neutron scattering data were taken at the NIST 
Center for Neutron Research. The integrated intensity 
of the (1, 1, 3/2) magnetic Bragg peak (Figure 3c) was 
measured using the BT7 spectrometer. The incident neu- 
tron energy was fixed at 13.46 meV and the horizon- 
tal collimation sequence was open-open-sample-40'-open. 
The (002) reflection of a pyrolytic graphite (PG) crys- 
tal was used as an analyzer. For the two-axis (energy- 
integrating) measurement of the instantaneous spin cor- 
relation length at T=70 K, we used the BT9 spectrometer 
with the sample aligned in the (HKO) scattering zone. 
The incident neutron energy was 35 meV, with a col- 
limation sequence of 40'-24'-sample-20'. In both mea- 
surements, a PG filter was placed in the incident beam 
to remove higher order neutrons. The inelastic neutron 
scattering data shown in Figure 4 were taken using the 
SPINS NG5 triple-axis spectrometer. The final neutron 
energy was fixed at 5 meV and the horizontal collima- 
tion sequence was guide — 80' — S — AO' (radial) — open. 
A liquid-nitrogen-cooled beryllium filter was placed in 
the scattered beam to remove higher order neutrons. 
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